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We study the monopole oscillation in the bose-fermi mixed condensed system by performing the 
time-dependent Gross-Pitaevskii and Vlasov equations. We study the resonant oscillation where the 
intrinsic frequencies of boson and fermion oscillations are same. 
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Over the last few years there has been significant progress in the production of ultracold gases, which 
realize the Bose-Einstein condensates (BEC) 

El HQ 

, degenerate atomic fermi gases |4| and bose-fermi 
mixtures 5]. These systems offer great promise for studies of new, interesting-driven phenomena. 
Collective oscillations are some of the most prominent phenomena common to a variety of many-body 
system. Theoretical studies on the collective motion in bose-fermi mixtures were performed in the normal 



3. 



phase with the „nle approach H and the random-phase approximation (RPA) 
In Ref. y| we constructed a transport mode for the bose-fermi mixing system by combining the time- 
dependent Gross-Pitaevskii(GP) equation and the Vlasov equation, and studied the monopole oscillation 
by solving the time-evolution of the system directly. Then we found a rapid damping in the fermion 
oscillation and a blurring of the beat in the boson oscillation. We considered the system where the boson 
number Nb is much larger than the fermion number Nf (Nb 3> Nf), which optimizes the overlap of bosons 
and fermions. In such systems the boson intrinsic frequency is -\/5 which has been shown in the system 
including the condensed bosons only |2|. On the other hand the fermion intrinsic frequency becomes 
larger as the attractive boson-fermion coupling increases. This is suggestive that the fermion intrinsic 
frequency becomes equivalent to the boson one for a certain value of the boson-fermion coupling. 
In this work we choose the value of the boson-fermion coupling where the boson and fermion intrinsic 
frequencies are same. Then we investigate resonant behavior of the monopole oscillation by performing 
the numerical simulation. 

Here we briefly explain our formalism. First we define the Hamiltonian for boson-fermion coexistent 
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system as follows. 

H = H B +H F + H BF (1) 



with 



Hb = J d\{-^\x)V 2 ^x)+ l -x 2 ^(x)^x) + 9 -f(^4>) 2 }, (2) 



h 2 i 

H F = I d 3 x{-- ^ty 2 ^ + -m f uj 2 f x 2 ^ip}, (3) 

H BF = h BF f d 3 x{^4npH}, (4) 



where <j) and ip are boson and fermion fields, respectively. Fermion mass mf and the trapped frequency 
ujf are normalized with the boson mass Mb and the boson trapped frequency Qb-, respectively. The 
coordinates are normalized by £s — (h/MsQs) 1 / 2 . The coupling constants g B and Kbf are given as 

9 b = ^waBB^B 1 ) ( 5 ) 
h BF = 2Tra BF (, B 1 (l + m~ f 1 ), (6) 

where a BB and a BF are the scattering lengths between two bosons and between the boson and the 
fermion, respectively. 

The total wave function |$(t) > has N c condensed bosons, whose wave function </> c is defined as 

4> c (x,t)=<&\4b(x,t)\$>, (7) 

where t is the time coordinate normalized with Qg 1 . In this work the wave function cf> c is expanded with 
the harmonic oscillator wave function u n (x) as 

N bass 

<f> c (x,r)= A n e^u n (bx 2 )e-^ x \ (8) 

n=0 

where Ni, ase + 1 is the number of the harmonic oscillator bases. We define the Lagrangian with the 
collective coordinates as 

L(A n , 6 n , 6, v) =< Hr)\{i^- - H}\<P(t) > . (9) 

Instead of solving the time-dependent GP equation directly, b, v as time-dependent vari- 

ables and solve the Euler-Lagrange equations with respect to these variables. 

The many fermion system can be described with the Thomas-Fermi approximation, which is given in the 
limit h — > 0. Here we define the phase-space distribution function as 

f(x,p,r) = J d 3 z < <f>\ip(x + ^z,t)^(x- ^z,t)|$ > e~ i P z . (10) 



In this classical limit this phase-space distribution function satisfies the following Vlasov equation 



y, 



^f(x,p;T) = (|- + -V, - [V x U F (x)][V P ] \ f(x,p;r) = 0. (11 i 
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with 



Uf(x) = ~ 



mfujjx 2 + HbfPb(x), 



(12) 



In the actual calculation we solve the above Vlasov equation 1)11(1 with the test particle method yj. 
Thus we can get the time-evolution of the condensed bosons and the fermions by solving the above 
time-dependent GP equation and the Vlasov equation. 

We calculate the monopole oscillation of the system 39 K — 40 K; the number of bosons ( 39 K) and fermions 
( 40 K) are taken to be 100,000 and 1,000, respectively. The trapped frequencies are taken to be = 100 
(Hz) and w/ = 1. The boson-boson interaction parameters is set to qb = 1.34x 10~ 2 which corresponds to 
o-bb = 4.22(nm) 6J; in addition we vary the boson-fermion interaction parameter Hbf which is measured 
in the unit of ho = 7.82 x 10~ 3 corresponding to asF = 2.51(nm) 
Using the root-mean-square radius (RMSR) R, we here define 



where Rq is the RMSR of the ground state, and r is the time variable normalized by fig . In Fig.H we 
show the time-dependence of this observable for bosons (Axb) and fermions (Axf) with hsF = — ^o- 
In this calculation the initial condition is taken to be Axb = and Axp = 0.1 at r = 0. We see a fast 
damping in the fermion oscillation as shown in Ref. 8]. However the the boson oscillation grows longer 
and, its amplitude becomes larger with Iibf — —ho than when the boson-fermion coupling is repulsive, 
h B F > 0. 

When varying the boson-fermion coupling, the situation becomes different. In Fig. we plot the results 
with Iibf = — 2h,Q. The oscillation behaviors are quite different from those in Fig. We can see the 
extension and shrink of the fermion density. 

In order to study this phenomenon, further, we examine the spectrum which is obtained by the Fourier 
transformation 



where Ac is a normalization factor. In Fig.[3]we give the spectra for boson (a) and fermion (b) oscillations 
obtained by integrating over the regions < r < 200. We can see that these two spectra have peaks at 
the same frequencies. 

The oscillation process is explained as follows. At first the fermion oscillation starts and trigger the 
boson oscillation. The boson oscillation grows while the fermion oscillation very rapidly damps. Once 
the fermion oscillation is stopped and grows up, again. The boson oscillation has a large amplitude 
at that time and scatter many fermions to the outside region. The boson density is distributed in a 
smaller region than the fermion density. Outside of the boson density the fermions feel only the harmonic 
oscillator potential, which is weaker than attractive force inside the boson density region. In addition the 



Ax(t) = R{t)/R - 1, 



(13) 




(14) 
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scattered fermions lose the kinetic energy and remain there rather long time. So fcrmions move from the 
inside region to the outside region in the expansion time 30 > t > 150. 

Later fermions far from the boson region return to the central region, but some fermions moving out 
from the boson region continue to populate outside the boson region because a test particle in the 
harmonic oscillator potential moves on closed orbit, which does not cross in the boson region. The 
system approaches to being equilibrium arounr r rs 200. Hence the averaged radius of fermion density 
remains to be 10 % larger than that at the ground state 

In Fig. 01 we show the intrinsic frequencies of the boson, 0J b M , (crosses) and fermion, (*> M , (full squares) 
oscillations as functions of the boson- fermion coupling hsF normalized with ho- For reference we also 
plot the fermion intrinsic frequencies with the open circles which are connected by the solid line, and 
u b M = \/5, which has been shown in the system only including the condensed bosons Q , with the dashed 
line. The boson intrinsic mode is almost independent of fermions because Nf, 3> Nf. 
When the boson-fermion coupling is Iibf ^ ~~ ho the intrinsic frequencies of the fermionic oscillation are 
almost equal to the frequency under the boson motion frozen, and their values are estimated to be Q 

w f « 2(w + & BF ^f ). (15) 

The frequencies are determined by the boson density distribution, and these fermion oscillations are 
volume vibrations. 

When the boson-fermion coupling is — 3ho < Hbf < —1.5ho, however, ujL, ~ 0J b M , and the resonant 
oscillation occur. This oscillation behavior is similar in the interaction region, though Axp reaches its 
maximum earlier as the interaction becomes more strongly attractive. 

When the boson-fermion coupling becomes larger, hsF ^ — 4/io, the fermion intrinsic frequency becomes 
u> M w 2.1, which is close to that only in the harmonic oscillator potential, u> M = 2. In strong attractive 
interactions the fermion motion loses collectivity, and RMSR of fermion density is affected only by the 
motion of fermions outside the boson density region. Namely the fermion oscillation is surmised to become 
a surface vibration. 

In summary we study the collective monopole oscillation in the bose-fcrmi mixed condensed system when 
fiBF = —2ho. In this coupling the intrinsic frequency of the fermion oscillation, u> M , is almost equal to 
with the intrinsic frequency of the boson oscillation, uj h Ml and we observed resonant behavior. In this 
resonant oscillation the amplitude of the boson oscillation becomes quite large even if the initial amplitude 
is taken to be zero. Furthermore the fermion density is largely extended because some fermions must 
move out from the boson density region to achieve the mixed system equilibrium. 

Furthermore we find that the fermion oscillation becomes a surface vibration when the boson-fermion 
coupling is more strongly attractive hsF ^$ — 4/joj than that in the resonant region, — 3/io < hsF 
— l.bho- While fermion oscillation exhibits the volume vibration when the coupling is attractive weaker 
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or repulsive, hsF ^ — l/io- I n future we would like to examine this consideration more detailed. 
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Figure 1: Time evolution of Aa; for boson (a) and fermion (b) at the initial condition that Sb = 1.0 and s/ = 1.1 
with the boson- fermion coupling hsF = — ho- 
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Figure 2: Time evolution of Aa; for boson (a) and fermion (b) with the boson-fermion coupling Hbf — —2ho- The 
other conditions are same as those in Fig. Q 
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Figure 3: Spectra of the boson (a) and fermion (b) oscillations shown in Fig. [5] 
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Figure 4: Intrinsic frequencies of boson and fermion oscillations versus the boson-fermion coupling. Crosses and 
full squares represent the intrinsic frequencies of boson and fermion oscillations, respectively. The open circles, 
which are connected with the solid line, indicate the fermion intrinsic frequency with the boson motion frozen. 
The dashed line denotes u)m = y/E. 



